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A first step towards a global description of a model that cannot be solved exactly is a mean- field theory (MFT). 
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A general expansion scheme based on the concept of linked cluster expansion from the theory of 
classical spin systems is constructed for models of interacting electrons. It is shown that with a suit- 
able variational formulation of mean-field theories at weak (Hartree-Fock) and strong (Hubbard-III) 
coupling the expansion represents a universal and comprehensive tool for systematic improvements 
, of static mean-field theories. As an example of the general formalism we investigate in detail an 

■ analytically tractable series of ring diagrams that correctly capture dynamical fluctuations at weak 

coupling. We introduce renormalizations of the diagrammatic expansion at various levels and show 
how the resultant theories are related to other approximations of similar origin. We demonstrate 
that only fully self-consistent approximations produce global and thermodynamically consistent ex- 
tensions of static mean field theories. A fully self-consistent theory for the ring diagrams is reached 
£SJ , by summing the so-called noncrossing diagrams. 
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It usually shares many features with an exact solution, offers a global phase diagram, and enables to investigate 
the existence of various solutions that exhibit long-range order. A comprehensive MFT represents a conserving and 
thermodynamically consistent approximation, i.e. there exists an explicit free energy functional applicable in the entire 
• range of input parameters. 

In the theory of interacting electrons it was for a long time only the weak-coupling Hartree-Fock solution that fulfilled 
the usual requirements imposed on a MFT in classical statistical mechanics of lattice spin systems . Attempts to 
construct a strong-coupling MFT around the atomic limit with analogous comprehensive properties failed in most cases 
to produce a thermodynamically consistent MFT. In the case of the Hubbard model mean- field theories reproducing 
the atomic limit have a structure of the Hubbard-III solution || and are commonly called alloy-analogy approximations 
, |^|- These approximations suffer from the deficiency that they do not represent a conserving approximation in the 
sense of Baym |5j. Only lately in connection with the limit of high spatial dimensions, d — > 00, a thermodynamically 
consistent extension of the Hubbard-III solution with a free-energy functional was proposed p] . 

However, having a thermodynamically consistent MFT alone is insufficient to make the mean-field results reliable. 
Unless we have an expansion scheme around the MFT with which we can study the stability of the mean-field solution, 
we cannot really judge whether and in which limit the mean field produces physically relevant results. Hence it is part 
of the construction of a comprehensive MFT to find an appropriate perturbation theory the first term of which is just 
the mean-field grand potential. It is fairly known that the Hartree-Fock MFT is a first, self-consistent approximation 
of the weak-coupling perturbation expansion. It is, however, much more difficult to find an expansion scheme the first 
step of which would be a MFT of the Hubbard-III type, since this MFT cannot be mapped onto a Fermi gas. We 
cannot use a usual expansion relying on Gaussian integration and the Wick theorem. 

We are hence facing the problem how to connect a general MFT with an appropriate perturbation theory. First 
attempts to abandon the weak-coupling expansion in the Hubbard model were made by Hubbard |?J who formulated 
an unrenormalized expansion in the hopping amplitude within which he rederived an approximation now called 
Hubbard-I. Recently Metzner [|| using a formulation in d = 00 extended the Hubbard approach and showed a way 
how renormalizations from the linked cluster expansion |^] can be introduced. However, it has proved cumbersome 
to deal with renormalizations effectively because of the quantum character of the model. This expansion does not 
efficiently go beyond the Hubbard-III solution and in this form it seems inadequate to improve strong-coupling mean- 
field theories. 

It is the aim of this paper to formulate an expansion scheme suitable to encompass general mean-field theories, the 
Hartree-Fock at weak coupling as well as the Hubbard-III at strong coupling, which enables to introduce consistently 
renormalizations leading to conserving approximations in the Baym sense. Importance of a systematic expansion 
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around a comprehensive MFT was recently acknowledged in M, where a formal perturbation expansion resulted 
from a specific variational formulation of mean-field theories being either an upper bound (Hartree-Fock) or a lower 
bound (extension of Hubbard-III from Q). In the following we formulate this expansion explicitly, show how to 
generate various approximations, how to introduce renormalizations, and study quantitatively the effects of dynamical 
fluctuations on the results of static mean-field theories. We use the linked cluster expansion (LCE) from the classical 
statistical theory of liquids and lattice spin systems reviewed in ^ . The LCE is applied here in quite a different manner 
than in ||, since we cannot expand in the hopping amplitude. We formulate the mean-field theories in such a way 
that we gain a universal expansion scheme for the weak- and strong-coupling theories. Within the derived scheme we 
compare different approximations having closed form. We further discuss the role and importance of self-consistencies 
in the light of the recent claim that a partial self-consistency, called iterated perturbation theory, may be superior to 
the fully self-consistent approach at intermediate and strong coupling [|l0| . We conclude that only fully self-consistent 
theories represent conserving approximations. As the main result we then derive a fully self-consistent version of the 
ring-diagram approximation by summing the so-called "noncrossing" diagrams for the grand potential and thereby 
win a global, thermodynamically consistent theory containing dynamical fluctuations that is equally applicable to 
weak as well to strong coupling MFT's. In this way we remove the doubts that renormalized RPA or more advanced 
self-consistent approximations may not lead to conserving approximations JTl| , ^2| . 

The paper is organized as follows. In Sec. II we reformulate the variational approach to the construction of the 
comprehensive MFT from Q so that the same expansion scheme can be used for the sample mean-field theories of 
interacting electrons: Hartree-Fock and the thermodynamically consistent extension of Hubbard-III from ||. Sec. Ill 
summarizes the formalism of the LCE needed to generate the appropriate expansion for the grand potential of itinerant 
models with local electron-electron interaction. In Sec. IV we choose a special class of diagrams, ring diagrams, leading 
to an approximation in closed form and being exact up to second order in the interaction strength. In Sees. V and 
VI we derive one-particle properties from the variational grand potential and present numerical results obtained from 
various simplifications of the ring-diagram approximation. Here the position of the iterated perturbation theory from 
pof in the present scheme is thoroughly discussed. In Sec. VII a controlled way to reach full self-consistency and 
thereby a thermodynamically consistent and conserving theory of the ring-diagram approximation is presented. 



A MFT is controllable only if it enables a systematic expansion towards the full solution. Not all constructions 
of mean-field like approximations are suitable for further corrections via perturbation expansions. E.g. heuristic 
derivations of the Gutzwiller extension of the Hartree-Fock MFT Jll| or of a strong-coupling MFT J2| do not seem 
to fit into a systematic expansion scheme starting with the MFT. On the other hand, a variational formulation of 
mean- field theories from jl| offers a consistent way to connect a MFT with a perturbation expansion. We present in 
this section an analogous variational formulation of mean-field theories of itinerant electrons that enables to formulate 
a universal expansion scheme for both weak and strong coupling. Explicitly we will refer to the weak-coupling 
(Hartree-Fock) and strong-coupling (Ref. Q) theories of the Hubbard model. 

For an expansion around a MFT we must decompose the total Hamiltonian into its unperturbed part (leading to 
the desired MFT) and a perturbation 



Our aim is to find suitable constituents of the decomposition (|lj) for the weak and strong coupling mean-field theories 
leading to conceptually the same perturbation expansion. 

We start with the simpler case of the Hartree-Fock approximation. For the sake of simplicity we restrict ourselves to 
spatially homogeneous solutions only. The variational parameters in the Hartree-Fock theory are the particle densities 
n a and we can write 



II. VARIATIONAL FORMULATION OF MEAN-FIELD THEORIES 



H = H + AH. 



(1) 



H — Hhf + AHhf 



(2a) 



with 




(2b) 




(2c) 
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where nj CT = c. 'c\a and M is the number of lattice sites. This decomposition leads to the standard perturbation 
expansion around the Hartree-Fock solution generated by the Hamiltonian Hhf- 

The situation in the case of the strong-coupling MFT is more complicated. The decomposition used in jl| leading 
to a lower bound does not have the structure of (0) . To obtain a structure similar to (||) we have to distinguish the 
static and dynamic degrees of freedom used in the strong-coupling MFT. We decorate the static degrees of freedom 
with an index s. We choose a decomposition of the Hubbard Hamiltonian into two Falicov-Kimball subhamiltonians 
with equal weig ht A = 1/2 0. We hence define 

H a = ( £ k + E a ) al a a ka - E^ a ^ n\_ a + ^ ^ ^K-a > ( 3a ) 

k i i 

where n\ a = cf^c^ and 

AH = Y,Ea (n| CT - n ia ) + - {^i-a - n|_ CT ) (3b) 

from which the desired decomposition of the Hubbard Hamiltonian directly follows: 

H = Y^ H ° + AH • ( 3c ) 

a 

It is easy to show that F — — lnTrexp{— /3J2a in d = oo generates the mean- field theory of ]jj with 
X a t = 0, A CT = 1/2. The variational parameters here are not n s ia (which are operators), but energies E a . It is also 
clear from (|3a|) that the dynamics of and decouple and Fmf = + Fj,. The stationarity conditions for such 
a MFT lead to 



( n icr)MF ~ ( n ia)MF • (A) 

We have two different contributions to the perturbation AH in (3b), one quadratic and one quartic in creation and 



annihilation operators. This makes perturbation expansion cumbersome. We can, however, drop the quadratic term 
without changing the physics of the full Hubbard Hamiltonian. The quadratic term only shifts the spin-dependent 
chemical potential [i a = [i + ah of the mobile (dynamic) electrons with the energy E a playing the role of the chemical 
potential for the local (static) electrons. These auxiliary local electrons completely decouple from the physical, 
itinerant electrons in the full solution. The energies E a in the physical sector then redefine the origins for the 
chemical potential /i and the magnetic field h. The physical results do not differ from those obtained from the original 
Hubbard Hamiltonian without energies E a . We can hence choose 

H = H MF + AH MF = J2 H ° + Ah mf (5a) 

and 

AH M F = — ^2n ia (r?i_CT - . (5b) 



i . a 



It is worth noting that the decomposition (^|)-(B) enabling a systematic perturbation expansion is not derivable for 
the MFT with Act = 1 and A at = — 1 as used in [| ,Q . To apply decomposition ([[]) consistently, no negative terms are 
allowed. If we choose X a t = 0, the resultant MFT does not reproduce the atomic limit exactly. Although the atomic 
limit is not covered by (]^), due to the halving of the interaction strength in the subhamiltonians H a , (JHJ) still contains 
the split-band limit. Important features of (ra) are that the mean-field free energy is an exact lower bound to the free 
energy of the Hubbard model and we can expand consistently around the MFT. Moreover, it was shown earlier in JlJ 
that the choice \ a — 1 and X a t — —1 is equivalent at T — to X a — 1/2, X at — with doubled interaction strength U. 

Both the mean-field theories for weak and strong coupling now have formally identical structure and the corre- 
sponding perturbation expansion can be formulated in the same manner. We use the LCE well elaborated in the 
classical statistical mechanics B as an appropriate tool for our purposes. 



3 



III. GENERAL FORMALISM OF THE LINKED CLUSTER EXPANSION 



The aim of the LCE is to represent the free energy (grand potential) as an appropriate series generated by derivatives 
w.r.t. auxiliary external sources added to the unperturbed Hamiltonian. Contrary to the standard weak-coupling 
expansion, LCE does not rely on the Wick theorem and is hence appropriate for the construction of a systematic 
expansion around any MFT. The only restriction on applicability of the LCE is that the generating unperturbed 
functional (free energy of the unperturbed system with suitable external fields) has to be explicitly known. 

Our aim is to formulate a perturbation expansion around quantum mean-field theories using decompositions (j^) 
and (||). Although the unperturbed mean- field functional to (Q) is known in any dimension, we restrict the expan- 
sion to infinite dimensions, since only there the unperturbed functional to the decomposition (^) can explicitly be 
constructed. This restriction leads to a substantial simplification of the perturbation expansion, since the effective 
electron propagator in perturbation theory does not depend on momentum explicitly [ jD^ ]. 

The generating functional for a quantum LCE can be represented as 



5mfW, Ea] = jj InTrTexp \ - J dr 



H M f(t) - ^M t ) T^ia IT J 



(6) 



with T standing for the time ordering operator and Tr denoting the trace in Fock space. Functional (|6|) is related to 
the grand potential ^\} F = —/3~ 1 g^ I > F . The generating external sources were absorbed into the mean-field variational 
parameters, n a for the weak coupling MFT and E a and the chemical potential of the itinerant electrons \i a for 
the strong-coupling MFT. These parameters now become time dependent. They are coupled to the time-dependent 
particle densities which are the only operators appearing in the perturbation AH . 
The Hamiltonian in (H) can explicitly be written for the weak-coupling MFT as 



H H f(t) = ^ e k aj ccr (r)a kcr (T) + U ^ n -v( T ) Ow(t) ~ n a) + UNn^ni 
k,u 



and for the strong-coupling MFT as 

H M f{t) = ^ a lA T ) a ^( T ) + E E <y{ T ) (™«x(t) - n| CT (r)) 

k,, 



-E 

2 ^ 



(7a) 



(7b) 



The time dependency was added to the variational parameters only where necessary. In the limit d = oo we obtain 
for the Hartree-Fock MFT 







9HFi n <y} = J dl ~ 



Un ] n i 



E 



dep(e) [ln{d/dr + p a - Un_ a - e)} (r, r+) 



(8a) 



while for the strong-coupling MFT the generating functional reads 



g { $ F { l iv,E a }=Y J Jdrl J dep(e)[ln(d/di 



Ho 



E„ 



<0] (r,r+) 



+ [lnG CT ] (7 



) + K(r)<(r) + (1 - nl») ([ln( 7(T )] (r, r+) - /T 1 Ml - 



r(r))) 



In 7<r 



(T,T+)-p- l ]nn s _ a (T) 



(8b) 



where n a (r,T') := 6(t - T')n CT (r), E {t,t') := S(r - T r )E a (T), h„{t,t') := S(r - t')m<t(t) and 7 .(t,t / ) := 
G~ x (t,t') + S CT (r,r'). The single-particle Green function G CT (r, t') — — {Tc„(t)c^(t'))t and the corresponding 
self-energy E CT (r, r') were used. With these expressions we can derive the explicit LCE with the perturbation AH(t) 
if we use 



AH hf (t) 



Udni(T) fe T (r) ' 



(9a) 
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8 8 
8^{t) 8E~Jf) 



(9b) 



and represent the full potential in d = oo for the Hartree-Fock MFT as 



gHF = In 

and for the strong-coupling MFT as 



exp 



dr 



<5n;(r) <5n T (r) 



(10a) 



9M F 



111 



exp 




(10b) 



Note that the imaginary time is removed from the variational variables in the end of all calculations to restore the 
time homogeneity. 

The aim of the LCE is to define a direct expansion for the grand potential, i.e. to interchange the logarithm with 
the variational derivatives in formulas (|ic|). Since the derivatives are of second order, this is a tremendous task and 
leads to sophisticated rules for the construction of exact contributions in all orders of the perturbation AH [|| . We 
can, however, formally write down the full potential as 



9mf = exp 




5^{t) 8E- a {r) 



Sm_f{av, E a } . 



(11) 



In ( |lip the brackets with variational derivatives denote a linear operator which commutes with the derivatives and 
acts on the unperturbed functional g^F as 



6n a {r) 8E- a {r) 



s2 9mf{^' E t} , s 9mf{^' E ^} S 9mf{^^ e <^} 



,(0) 



Co) 



5^(r)8E- a (r) 



8E^(t) 



(12) 



Definition @ expresses a derivative w.r.t. a bond (nonlocal) variable (left-hand side) in terms of site (local) 

variables genuine to the unperturbed functional g^F (right-hand side). The power expansion of the exponential in 
( pd| ) together with the definition of the bracketed derivative ( [l2"| ) determines the unrenormalized LCE. To save the 
space we make explicit derivations for the strong-coupling MFT only. We hence drop the subscript "MF" in the 
forthcoming formulas. 

It is convenient as in other perturbation expansions to define a graphical representation for the LCE and to deal 
with diagrams by introducing vertices and internal lines. If we denote x"(t), a = ± with x~£(t) = /v(r) — Ha and 
x~ (r) = E a — E a (r) we can define an unrenormalized n- vertex function 



F (n) ai—a n /_ 

1 (0) CTl --<7„ V 1 ' 



8 n gW{y a ,E a } 
Sx^iri) ■ ■ ■ 8xZ{r n ) 



(13) 



and an unrenormalized internal line (edge) (r±, t-i) connecting the vertices. The "propagator" is hence a matrix 

in the internal degrees of freedom of the vertices, a, a and r. We can assign a graphical representation to both the 
vertices and the edges. We denote an n- vertex as 




(14a) 



where 1 = (ai, <7i, ti), ... . In this version of the LCE the electron-electron interaction plays the role of edges 
connecting the vertices. We assign a dashed line to an edge, i.e. 
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* — v%%(n,T5) := -6( Tl - T 2 )6 au - a3 6 ai ,- a3 . (14b) 
2 Z 

With this graphical representation we can try to sum various subclasses of diagrams either with or without renor- 
malizations |9|. The renormalized quantities are standardly denoted with full circles and bold lines, respectively. Note 
that in the present formulation of the LCE the vertices contain the one-electron propagators (Green functions) while 
the edges represent the electron-electron correlation. 

To obtain a consistent LCE we must treat all the variational functions and parameters that do not serve as generating 
sources in the unperturbed, mean-field functional as external variables independent of the generating fields. The mean- 
field stationarity equations, leading to an approximate grand potential of the problem, are imposed only in the end of 
the calculations after a chosen class of diagrams has been summed up and the time-dependence from the generating 
fields removed. It is hence necessary to have an expression for the grand potential to be able to expand consistently 
around a MFT. 
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IV. RING DIAGRAMS 



The LCE is of practical use only if we are able to pick up a class of diagrams that can explicitly be summed and 
then numerically evaluated. Due to the internal degrees of freedom (Matsubara frequencies) the quantum version of 
the LCE is much more complicated than its classical counterpart. We hence cannot expect to go that far in summing 
various classes of diagrams as in the case of classical fluids or spin systems. In the quantum case we are practically 
restricted to ^-derivable approximations (l6| where higher order vertices (n > 2) are neglected. We therefore resort 
to LCE diagrams with renormalized 1- and 2-vertices only, i.e. to vertices that are at most second derivatives of 
the renormalized functional g. An approximation for the functional g( N * containing maximally TV-vertices can be 
represented in the Baym form M 



exp 



' N 

5>[e 



9 {0) {x} 



x=0 



N 



(15) 



where T n is a renormalized n- vertex, C n is a renormalized n-edge and V x is a derivative in appropriate generating 
sources x. Symbols with hats stand for matrices. Here the trace tr is taken over all the internal degrees of freedom. 
The use of generally renormalized quantities is crucial here, since only then the functions F„ and C n can be treated 
as independent variational variables and the resultant approximation can lead to a thermodynamically consistent 
and conserving approximation. The matrix variables C n and T n are determined from stationarity equations, i.e. 

If we restrict the theory to N — 2 we are left with two types of diagrams contributing to <&. The former is the 
effective-field contribution involving only IV The functional $W then takes the form 



dr 



6^(t) 5E_ a { T y 



(16) 



Since we expand around a mean field where E a are variational parameters, the derivative Sg^ /8E- a = for arbitrary 
r G [0,/3]. The "effective-field" approximation around a mean-field must hence vanish. This is a consequence of 
stationarity of the mean-field approximation w.r.t. all variables. 

A further step and the first nontrivial correction to a MFT within the LCE are the ring diagrams consisting of single 
loops with only 2-vertices F2. It is easy to find a diagrammatic representation for the ring diagrams contributing to 
the functional $: 



G 



To find a mathematical expression for the sum of ring diagrams we must first explicitly write down a matrix 
representation for both the renormalized vertices (full circles) as well as for the bare edges (dashed lines). They 
are 4x4 matrices in the case of the strong-coupling MFT and 2x2 matrices in the weak-coupling MFT. In the 
strong-coupling MFT we denote T"*%* = S 2 g^ /6x"l5x"% and represent the 2-vertex as a matrix 



r 2 = 



r++ r+- 

aa' aa' 
aa' aa' 



(17) 



where the single matrix elements are 2x2 matrices. The bare edge v can be represented as another 4x4 matrix 



v= —5{t-t) 



ax 
a x 



(18) 



where ax = 



1 

1 



is the first Pauli matrix. Using this matrix representation we can write down the sum of ring 



diagrams for the grand potential as 



$ (2) [P] = --trm(i-uf; 



(19) 



The prefactor 1/2 compensates the doubling of contributions to the grand potential due to the chosen matrix repre- 
sentation. Inserting (|l9|) in ( |l5| ) we obtain 



,(2) 



cxp ^ -tr 



g 2 v 2 



3 (0) W 



- -tr In (l - vf 2 ) - -tr 
x =o 2 V V 2 



C 2 T 2 



(20) 



where C 2 and V x are also 4x4 matrices. Approximation (g(J) was used by Horwitz and Callen |l7| and Bloch and 
Langer jl8j to improve the Weiss MFT for the Ising model. Although solvable at the classical level, it is not viable at 
the quantum level where the generating sources are time dependent. The exponential in (20) leads to a new functional 
integration being as complicated as in the original problem. We hence have to approximate the exponential. The 
simplest way to do that is to use the power expansion in C 2 and to terminate it appropriately. Since the most important 
dynamical fluctuations come from second order perturbation theory, we resort to an approximation reproducing the 
exact result asymptotically at least to U 2 . A more complete treatment will be presented in Sec. VII. 

To reproduce second order of perturbation theory it is sufficient to take into account only the linear term in the 
operator C 2 . Then (EOT) reduces to 



(2) = g(0) + 1 tr 



cw 



9 {0) {x} 



- -trln (l-v? 2 ) - -tr 
x=o 2 V J 2 



(21) 



Such an approximation is already explicitly solvable and leads to an analytic expression for the grand potential 
fi( 2 ) = —j3~ 1 g( 2 \ We use the matrix representations for T 2 and v , (|l7j ) and ([l8|), and evaluate the second derivative 
of the unperturbed functional g 1 - ^ . It immediately follows from (|8b|) that the derivatives with different spin indices 
vanish and the matrix of the second derivative in (El]) substantially simplifies especially in the frequency representation. 
Stationarity of the functional w.r .t. C 2 and T 2 yields explicit expressions for these matrix functions. We exclude 
them from the potential g^ and consider the paramagnetic phase exclusively, so that the quantities loose their spin 
dependencies. We then obtain for the physical grand potential 



^oS-V = -En s - /r 1 e "" 0+ \ J <M e ) ln ^ + 1* ~ E ~ S ( ZCJ «) - £ 1 + ln G (^«) 



+ (1 - n s ) [ln 7 (^„) - /T 1 ln (1 - n s )] + 



In 



U 



/rMnn" 



E e^ 0+ ln[l-ac/{(l-^)(G( )G( ))^ m )+n s (G^G^)^ m )} 

1 q=±1 m= — oo 



(22) 



Here w„ = (2n+ l)n(3 1 , v m — 2mnfl 1 . We further used the abbreviations j(z) = G(z) 1 + (z) — 7(2) 1 , 

GM( Z ) = [ 7 (z)_ vy 1 and 
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oo 



(23) 



with a = 0, U. The mean-field parameters E and n s , and the complex functions H(iu) n ) and G(iu> n ) at fixed \i are 
treated in the approximate grand potential variationally in the same way as in the unperturbed theory with 
The new grand potential Q ^ defines then a thermodynamically consistent theory. 

An analogous expression can be derived for the expansion around the weak-coupling MFT. Explicitly we obtain 



2Af ~ 0. P 2^ 



1 

4£ 



-ln 7 ^„)} + ^ J2 E e " m ° + 



dep(e) In [iw„ + ji — Un — T,(iu) n ) — e] + In G(ito n ) 
\u 1-aU (G i0) G (0) ){iis m ) 



a— ±1 m- 



(24) 



where the function G^ has the same meaning as in (22 ). Grand potentials ( p2[ ) and (]24|) hold only for the paramagnetic 
phase. Grand potentials allowing for a broken symmetry can be derived without conceptual changes. 

Having explicit expressions for the contributions due to the ring diagrams we can compare the resulting th eory 
with other similar approximate schemes, especially those studied extensively by Menge and Miiller-Hartmann [p.9| , 
iterated perturbation theory (IPT) introduced recently in connection with the metal-insulator transition in 

the d = oo Hubbard model, and with the extended Edwards-Hertz approximation |22] , |23]| . This will be done in Sees. 
VI and VII. Before we do this we should realize that the presented sum of diagrams within the LCE scheme does not 
rcnormalize the one-electron Green function but rather the interaction strength. The ring-diagram approximation 
(p2|), ( p4[ ) represents a sum of a geometric series of 2-particle bubbles renormalizing the interaction, i.e. diagrams 
contributing to the polarization of the vacuum. Because of the locality of the electron-electron interaction in the 



Hubbard model only even powers of U contribute. Grand potentials (23), (24) correspond to an "effective- interaction" 
(bubble-chain) approximation studied at the level of self-energy in [19). However, the derived sum of ring diagrams 
does not represent a fully self-consistent approximation, since the local electron propagators in the expansion are 
unrenormalized, mean-field ones, G(°\G( u K Although these propagators contain the self-energy S and are not free, 
the dependence of G^ on S is caused by the reduction of a d = oo theory to a d — (single-site) problem. Such 
a partial self-consistency is of topological rather than of dynamical origin. Moreover, neither (^2|) nor (|24|) has the 
Baym form (f| where the ^-functional depends only on the full propagator G. We hence have no guarantee that they 
are conserving approximations. 



V. ONE-PARTICLE PROPERTIES 

Grand potential ( p2| ) leads to a numerically solvable theory. In the present form it is only suitable for the para- 
magnetic phase and we use it to evaluate explicitly the one-electron properties. In order to simplify the numerical 
procedure we restrict our considerations to zero temperature. 

The defining equations for the variational variables are derived from stationarity of f2*- 2 ) . Vanishing of the derivatives 
w.r.t. the static variables n s and E leads to 



OO 

E 

711 — — OO 



Q = ±l 



C a (iv m ) [(G*°>G<°>) (ivm) - (g^G^) {iv m ) 



(25a) 



OO 

OO „ 

p-i J- e --o + y dep(e)[iu n +fi-E-E(ioj n )-e]- 1 



(25b) 



TL— — OO 



where 



f = -r l E 



i-|g(°)^„) 



(25c) 



and 
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C a (z) 



l-aUi(l 



n s )^G^G^)(z) 



i s (g^g^ 



(*)} 



(25d) 



The term in (25a) including a sum over bosonic frequencies is a correction to the unperturbed static MFT due to the 
ring diagrams. Variations of fl^ w.r.t. H(iuj n ) and G(iuj n ) lead to corrected mean-field equations for these functions: 



oo 

G(iuj n ) = J dep(e) [iuj n + fx - E - £(««;„) - e] _1 



(26a) 



and 



G(iu n )= (f - n s )G {a \iLu n ) +n s G (u \iuj n ) 

(1 - n s )G^(tLO n ) 2 a E ei% ° +G(0) (^n + ^m) [Ca(^m) + C a (-iu m )] 



a— ±1 m— — oo 



4/3 

-— n s G^(^n) 2 E a E e^ 0+ G( c/ )^„ + i ^ m )[C ct (^ m ) + G Q (^ m )] . 

a— ±1 rn— — oo 



(26b) 



To facilitate the numerical evaluation of the derived equations we convert the sums over Matsubara frequencies to 
integrals over real frequencies. However, this is possible only if the functions under consideration are analytic in the 
upper and lower complex half-planes. Apart from the one-particle Green function, the analytic structure of which is 
clear, there are two other complex functions (GwGW) (z) and C a (z). The former function is analytic with a cut 
along the real axis as can be seen from the representation 



X (a \z) ■= (G {a) G (a) ) (z) 



-— I il.r 

7T 



G (q) (x + z) + G [a) [x - z)] ImG (Q) (x + i0 + ) 



(27) 



On the other hand the function C a {z) can have a pole in the complex plane if 

1 = aU{(l - n s )X°(z) + n s X^ u \z)\ . 



(28) 



Since X^{z) = X^{—z), the r.h.s. of (|28f ) is real for z = iy, i.e. along the imaginary axis. There may be a 
critical value U c at which the pole appears for the first time at z — and where the perturbation expansion breaks 
down. We must hence distinguish two regimes of the interaction strength. If U < U c (weak coupling) we have no 
polar contributions to the spectral representations of the Matsubara sums. At strong coupling, U > U c , the pole in 
C a (z) influences the integrals over real frequencies. The existence of a pole in a two-particle Green function indicates 
an instability of perturbation theory and is connected with a metal-insulator transition. We postpone a detailed 
investigation of the actual role of such a pole to a separate publication and stick in this paper excl usive ly to the 
weak-coupling regime U < U c . In this case it is easy to show that the sum over bosonic frequencies in (25a) vanishes 
due to analyticity of the integrand in the upper and lower complex half-planes . Th erefore the mean-field equations for 
n s and E are not changed by the contributions from the ring diagrams. Eq. ( [26b| ) can be rewritten at weak coupling 
for an arbitrary complex energy z as 



U 



G(z) = (1 - n s )G^(z) \ 1 + !±G (0 \z)I^(z) \ + n s G^(z) I 1 + ^G^(z)I^(z) 



where we used 

/<■>(*)= E 



a=±l 



dx 



C a {z - ir)IniG (a) (i + i0 + ) - G {a) {z + x)ImC a (-x + i0 + ) 



(29a) 



(29b) 



The function C a (z) is defined through the functions X^ as 



C a (z) = {l - aU [(1 - n s )X^(z) + n s X^ u \z)\ } 



(30) 
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Equation (29a) can be turned into an equation for the self-energy if we use the definition of G^ given below 
eq. (H|). We then obtain 



E(z) 



U 



2 1 + G{z) [E(z)-{7/2] 
U l-n s 



2 1 + G(z)T,(z) 
Eq. (|3l|), completed with the mean-field equations 

o 



I<°>(*). 



dxlm In 



and 



I l + G(z)S(z) {u) 
1 + G(z) [E(z) - U/2] K ' 



U 



G (0 \x + iO+) 



1 f 

-— / dxlmG(x 

7T J 



(31) 



(32a) 



(32b) 



is to be solved numerically for z = uj + i0 + . From (32b) we see that n s = n as is the case in the static MFT from 
pL Due to the electron- hole symmetry it is sufficient to solve the equations for ui < 0. To obtain the values for the 
self-energy and the one-electron Green function on the positive frequency axis we use the electron-hole transformation 
(n — > 1 — n, lo — » —ui) leading to: 

(33a) 



G(-u> + i0+) = -G(uj + i0+)* 



and 



E(-w + i0 + ) = - - E(w + i0+)* 



(33b) 



The chemical potential transforms in this expansion as /i — > {7/2 — /i and the parameter E 1 — » — _E. It is to be proved 
that (|33| ) really holds for a solution of (^l]), i.e. the approximate theory correctly exhibits the electron-hole symmetry. 
The proof is presented in the Appendix. 

The same steps can be made to derive equations for the self-energy of the ring diagrams summed around the 
Hartree-Fock grand potential (E4) . The resulting self-energy then is 



J HF 



(*) = 



{7/2 



l + G(z)E(z) 



J<°>(z). 



(34) 



We see that the self-energy from the strong-coupling MFT formally reduces to T^hf if we put n s = and E = Un. 
The electron-hole symmetry in the weak-coupling case takes the form 



G(-u> + i0+) = -G(ui + i0 + )* , 
Z hf (-lu + i0 + ) = -Y, HF {u + io+y 



(35a) 
(35b) 



VI. NUMERICAL RESULTS AND COMPARISON OF VARIOUS APPROXIMATE SCHEMES 

We now apply the derived equations in the weak-coupling limit to assess quantum fluctuations due to dynamical 
corrections to the static MFT. Only for the sake of simplicity we perform the numerical calculations on a Bethe lattice 
in d = oo, i.e. we use for the diagonal element of the one-electron Green function the representation 

-l 

(36) 

where D — 2 equals half the bandwidth. We also restrict the calculation to the half-filled band in the paramagnetic 
phase, i.e. n a = n% = 1/2 and /j,= U/4 (/i — E = 0). Furthermore we remain during all the calculations at T = Q. 

Since our aim is to investigate the resulting equations only in the weak-coupling limit, we separately analyze the U 2 
corrections to the static MFT before we take into account all the ring diagrams. In course of calculations we compare 
the results obtained from the LCE treated here with other approximations of similar origin. 



G^(z) 



G{z) 
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A. C/ 2 -corrections to the mean-field solutions 



Let us start with the simpler expansion around the Hartree-Fock MFT. Using ([m|) we obtain for the Z7 2 -contributions 
to the self-energy 



, 0) r-> } dx 



l + G(z)2J(z) u 1 + G(2;)2](z) 7 7r L 

— oo 

+G (0) (z + a;)ImX( )(a; + iO + )j . (37) 

The resulting self-energy consists of two different contributions, the algebraic prefactor and the integral Iq°\ The 
former is a consequence of the incomplete self-consistency of the defining equation (|37j). The Green function 
from the r.h.s. of (p7| ) is not identical with the full one-electron propagator G. Hence the variation of the grand 
potential (E2) w.r.t. G, leading to the equation for the self-energy, generates algebraic factors due to the variation of 



the propagator G^K The latter contribution to (37) is exactly the result of the iterated perturbation theory (IPT) 
and corresponds to a partially self-consistent, second-order Feynman diagram for the self-energy. It is important that 
both the contributions are complex and lead to a rather intricate analytic structure of the self-energy S(z). It is 
impossible to prove the desired Herglotz properties, i.e. X(z) = S(z*)* and ImE(z) oc — Imz for all complex energies 
z. Only at the Fermi level, u> — 0, we obtain from the electron- hole symmetry that E(0) = 0. 

Approximation ( |37| ) does not seem to cause troubles at weak coupling (U < 1.9), where only a slight modification 
of the non-interacting density of states (DOS) can be observed (Fig. 1). However, interactions stronger than U c ~ 1.9 
lead to a breakdown of the approximate equations, since a singularity appears in the prefactor [1 + G(z)E(z)] . 
Namely we observe that the derivative dG^\uj)/dui diverges at a critical frequency uj c , which leads to a jump in the 
propagator G^°\uj). The density of states then develops a cusp at the critical frequencies ±o> c (Fig. 2). In such a 
situation it is impossible to reach convergence and it is unclear if the equation ( |37| ) does have a solution at all. 

The algebraic prefactor is a consequence of the variation of the grand potential and has no diagrammatic background. 
It causes the breakdown of the iteration scheme. We may ask if we can do better by neglecting this prefactor. Such a 
choice corresponds to the diagrammatic expansion directly applied to the self-energy instead of to the grand potential. 
We obtain 

o 

S(z) = C/ 2 J ^ \ i X^\z~x)lmG w (x + tO + ) + G {0) (z + x)lmX ( ^(x + tO + )^ , (38) 

— oo 

which is the IPT self-energy. Analyzing (|3"s| ) we find that this equation can safely be iterated up to U w 4.7 (Fig. 
3) where the results of Georges and Krauth J2lj are reproduced. The DOS shows a three-peak structure where the 
central band narrows with increasing U. However, it seems impossible to reach the alleged metal-insulator transition 
at which the width of the middle band has to vanish. Before we reach this point, the iterations break down for similar 
reasons as in the solution of d37j ) . With increasing interaction strength a discontinuity in the real part of the Green 
function G?(°)(w) seems to develop together with a cusp in its imaginary part (Figs. 4a, b). The extremal points 
for both the functions approach each other (Fig. 5) and seem to merge at oj c < 0, i.e. away from the Fermi energy. 
Earlier a singularity in G^ was reported on the real axis only at the Fermi energy |^(|. According to our results the 
Green function G^ displays a discontinuity at uj c and it is again impossible to reach a solution via iterations. If we 
attempt to solve ( |38| ) for complex energies, the breakdown of the iteration scheme is shifted to stronger interactions 
and the critical frequency is pushed towards the Fermi energy. Although the singularity smooths in comparison with 
the real frequencies, we did not succeed to remove it completely from G^ and to reach a convergent solution at 
strong coupling. If there were a singularity at G^ even at complex frequencies, the Green function could not fulfil 
the desired Herglotz properties. However, if the Herglotz properties of G^(z) were violated, it would be impossible 
to continue analytically the quantities from the imaginary frequency axis to the real energies in a unique way. The 
approximation then would become unreliable with possible spurious poles (zeros) in the complex plane. It is hence 
of great importance now to prove or disprove the Herglotz properties of the IPT self-energy if we want to draw 
conclusions from this approximation especially at strong coupling. 

From the above results it may seem that the IPT for the self-energy be superior to the LCE for the grand potential 
because of the unphysical complex prefactor in (p7j). This prefactor appears due to the fact that the one-electron 
propagator G^ , used in the ring diagrams, deviates from the full propagator G and thus the expansion is not fully self- 
consistent. It is a consequence of the non-self-consistent character of the expansion that not all physical quantities 
have a diagrammatic expansion. Diagrammatically the derivative is not connected with a cut of the one-electron 
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propagator only, but a variation of G^ w.r.t. G appears. This variation leads to the lost of the diagrammatic control 
over physical quantities and causes violation of Ward identities. Undesired prefactors are inevitable companions 
of non- or only partly self-consistent theories. Although nonexistent at the level of self-energy within the IPT, the 
unphysical prefactors do appear in higher order Green functions and correlation functions (susceptibilities). Therefore 
we cannot rely on IPT susceptibilities pl[ . Since the existence of "spurious" complex prefactors is a ubiquitous feature 
of partly self-consistent theories, we may conclude that only fully self-consistent theories in the spirit of Baym offer a 
physically consistent global picture. 

We can also calculate the U 2 corrections to the strong-coupling MFT. There we meet the same problems with 
prefactors and a breakdown of the iteration scheme as in the case of the expansion around the Hartree-Fock MFT. 



B. Corrections to the mean-field solutions due to the ring diagrams 



One may still hope to improve the weak-coupling approximation and to avoid the breakdown of the iteration 
scheme by taking into account the ring diagrams. However, the theory with ring diagrams suffers from the same 
deficiency as the U 2 terms, namely the equations contain unphysical complex prefactors, and an analogous breakdown 
of the iteration procedure is to be expected. It turns out that for the ring diagrams a singularity in the prefactor 
[1 + G(z)T,(z)] shows up even for smaller values of the interaction, namely U « 1.66, where the approximation 
scheme breaks down. Above this value the iterations do not converge any longer. It indicates that the additive 
diagrams push down the "transition" from the weak to the strong coupling regimes. To confirm this we suppressed 
the prefactors and thus gained an extension of the self-energy from the IPT (|38| ) due to the ring diagrams. The result 
is plotted in Fig. 6. We can see that the self-energy has a narrow central peak with a tendency to build up two satellite 
bands. However, before the satellites can be fully formed a pole in the two-particle function C a {oj) (30) at uj = is 
reached at U ~ 1.9 and the theory breaks down. There is no way to go around the pole in the two-particle function 
C a {z) in this not fully self-consistent approximation. The pole causes the Herglotz properties to be manifestly broken. 
There is no convergent solution for U > 1.9. The only way to restore analytic properties of the local, one-electron 
Green function is to make the theory fully self-consistent. 



VII. SELF-CONSISTENCIES AND RENORMALIZATIONS IN THE LINKED CLUSTER EXPANSION 

We saw in the preceding section that non-self-consistent diagrammatic approximations lead to differences in physical 
quantities derived either diagrammatically or with the use of derivatives (variations) of the grand potential w.r.t. 
external sources. This is a general feature of not fully self-consistent approximations where some Ward identities are 
violated. This conclusion applies also to the approximations investigated in Sees. IV- VI. Although the corresponding 
grand potential has a form where all the unknown functions have to be determined from stationarity conditions, the 
internal electron propagators in the diagrams are G^ — (G _1 + S) _1 (or G^ u > = (G _1 + S — U /2) -1 ) instead of G as 
one expects in the Baym approach. These are dynamically free propagators, since the dependence on the self-energy 
does not stem from the dynamics but rather from a special topology of the d = oo lattice mapped onto a single site 
(impurity) embedded in a medium. E.g. if the diagrams were treated in d = 3, the function G^(z) had to be replaced 
by the bare (unperturbed) k-dependent propagator G^(k, z). There would be no self-consistency at all in hnite 
dimensional lattices. The origin of the partial self-consistency in the investigated approximate schemes lies entirely 
in the reduction of the topology of the d = oo model to a "self-consistent" single-site problem j24[] . The topological 
self-consistency in the local propagator G^\z) coherently copies the dynamical processes from the impurity to the 
medium which, on the other hand, provides the electrons for the impurity. 

In this light it is clear that e.g. the self-energy, obtained by varying the grand potential from the LCE, does not 
equal the self-energy obtained from a sum of the non-self-consistent self-energy diagrams. To get rid of such ambiguity 
problems we must abandon non-self-consistent theories and introduce true dynamical renormalizations into the LCE. 
We can try to introduce the self-consistency naively via a replacement of the unperturbed electron propagator G^ 
with the full propagator G. However, such an introduction of the self-consistency into the theory is dangerous, since 
it need not be diagrammatically controlled. This naive way of turning the theory self-consistent may lead to multiple 
counting of diagrams and to undesired unphysical and spurious structures. The only consistent way to introduce 
self-consistencies into diagrammatic expansions is to classify diagrams into particular classes which may be summed 
by iterations, i.e. by self-similar insertions. 

Possible renormalization schemes within the LCE were extensively discussed by Wortis 11. However, in the quantum 
case with internal time (frequency) degrees of freedom, we are much more restricted to produce a tractable approx- 
imation. As we already mentioned, even a quantum analogue of the Bloch-Langer approximation from the classical 
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Heisenberg model [Q is not viable in the quantum case. Due to the time dependence of the vertices the ordering of 
functional derivatives, i.e. the time ordering of the external legs at vertices, is important. Interchange of derivatives 
leads to "crossing" of interaction lines and to different quantitative results for diagrams with vertices of higher order 
than two. If we resort to a fixed time ordering of the vertex legs, i.e. to a kind of "noncrossing approximation" we 
obtain a self-consistent version of the ring-diagram approximation. 

To this end we must take into account noncrossing diagrams resulting from the expansion of the exponential 
cxp |^tr C2V^j | in (pp|). If we analyze the contributions from products of the operator G2V 2 we can classify them 

into crossing and noncrossing diagrams, according to whether the two-point function G2 is or is not crossed by an 
interaction line (Fig. 8), respectively. The topology of higher order diagrams becomes more and more intricate due 
to the crossed diagrams (Fig. 8d) . If we neglect them we can exactly sum up the contributions from the noncrossing 
diagrams via self-consistency. Insertions as in Fig. 8c can be summed by replacing G^ ' with G, where the self- 
energy due to the C-contributions (hatched bubbles) equals to a convolution G * C. The grand potential due to the 
self-consistent ring diagrams around the Hartree-Fock MFT can then be written as a generating Baym functional 

f °? 

a n=-oo [-oo 

-AT, a (iu n ) - e] + A£ ff (wj„)G ff (MJ„) + In [l + G a (iuj n ) (S,(iu„) - (G cr C cr }(iLJ n ))] } 
+/T 1 f] e w '" 0+ ^T a {iv m )C a {iv m ) + i In [l - U 2 T^iu m )T i{iv m )] j (39) 



m— — 00 



where we used the same notation for the convolution {G a C a )(iuj n ) as in (p3|), but with the sum over bosonic Matsubara 
frequencies. Here AS^ is a correction due to the G,j-induced contributions from the noncrossing diagrams and G a is 
its Legendre conjugate variable. If they vanish, (|39]) reduces to its non-self-consistent version (p4|). If we exclude the 
variables r CT , C a and G^, AE^ we can reduce (|39|) at the saddle point to 

00 ( 00 

\-Sl%™ 9) = -Un m -r'Y, E e ""° + \ I d( P^ ln ^« + ^ - Un- a - E ff (ia; n ) - e] 



a n — — 00 



1 00 

+G CT (iu; n )£ CT (iu; n )} + — £ e iv ™ 0+ In [l - U 2 (G r G r )(iv m )(G L G L )(iv m )] . (40) 

m— — 00 

Grand potential f^™ 9 -* is the fully self-consistent "bubble-chain" approximation investigated at the level of self-energy 
by Menge and Mullcr-Hartmann [ fL9| . 

Analogously to (^) we derive a self-consistent version of the ring-diagram grand potential for the LCE around the 
strong-coupling MFT: 



00 



a n— — oo 



+ (1 - n s _a) ln [l + G a (iuj n ) (T. a (iuj n ) - (G a C a )(iu n ))] + A^ a (iuj n )G a (iuj n ) 
+nf_ CT ln 1 + G a (ioj n ) ^£ CT (iw„) - ^ - (G CT G CT ) (iu%)^ } ~ E 

+/T 1 £ e w ™ 0+ |^r CT (i I / m )G CT (w m ) + iln [l - [7 2 r T (w m )r ; (M/ m )] 1 . 

m— — 00 ^ ct J 



(41) 



Excluding the variables r CT and C a at the saddle point we obtain defining equations for the Green functions G a and 
G a 

oc 

G a {iuj n ) = G a (iu n ) = / dep(e)- ■ — — — r — — — — r (42a) 

J iuj n +fi-E tr - Y, a (iu) n ) - AS CT (iuj n j - e 

—00 

and for the self-energies Yi a (iwn) and A£ CT (ta; n ) 
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AE CT ^„) = -f-J £ G a {iu n + iv m ) - TT2 fr^ } [Trrv- \ ' ( 42b ) 

U n s _ 
E a (iui n ) - AEviiUn) = — ,. — w — Av ,. — ; — 777^- • (42c) 
2 1 + G ff (jw„) {^(lUn) - AL a {iuj n ) - U/2) 

Grand potential Q^p 9 ^ has an attractive feature in that it contains the split-band limit of the Hubbard-III solution 
and is exact at weak coupling up to U 2 . In contrast to previously studied theories with analogous properties p2|,f23[, 



grand potential ( J4l| ) is derived in a fully controlled manner. There are no spurious prefactors in the defining equations 
for the self-energies calculated from ( p9| ) and (|4l]). There is also no artificial breakdown of the iteration scheme in 
the fully self-consistent theory. All the Ward identities are now fulfilled and there is no difference in deriving various 
physical quantities either from varying the grand potential or using the corresponding diagrammatic expansion for 
them, since the theory is thcrmodynamically consistent and conserving in the Baym sense. 

We can now analytically continue the functions E a , AE CT and G a to the whole complex plane and use spectral 
representations with real frequencies. The numerical solution is then obtained straightforwardly. The quantitative 
results in the paramagnetic phase (Fig. 8) obtained from ( [40|) at weak coupling (U < 3.5) agree with the results of the 
"bubble-chain" approximation from pj| . For stronger coupling a pole in the two-particle function C a (uj) interferes 
and the iteration scheme must be modified to preserve the Herglotz properties of the one-electron propagator |25| l . 

Self-consistent mean-field theory with dynamical fluctuations ( [42] ) exhibits similar behavior at weak coupling (U < 
1.5) (Fig. 9). The DOS at the Fermi energy begins to decrease, since the theory is not a Fermi liquid (the term 
proportional to U 3, violates the Luttinger theorem). Two satellite peaks start to develop and split off at about 
U c w 5.5 (U c — 4 for the static MFT) similarly to |2^]. Because of frequency convolutions there are no sharp 
(algebraic) band edges and it is numerically hard to determine when exactly the bands split off. We also had troubles 
to reach convergence for the frequencies lying within the expected gap. To reach a reliable solution in the insulating 
phase, the iteration scheme must be refined. We, however, did not meet any singularities indicating an instability 
of this dynamical mean-field theory. The two-particle function C a {ui) did not reach the pole ( p8| ) because of the 
metal-insulator transition and we did not observe any indication for a divergency in the perturbation theory. 

The self-energy from ([l2]) contains two parts, AE being the direct contribution of the ring diagrams and E — AE 
which is the mean-field self-energy with the renormalized one-electron propagator. To assess these two different 
self-energies we plotted their imaginary parts separately in Figs. 10a, b and compared them with the corresponding 
self-energies from the Hartree-Fock solution ([l0]) and the static MFT |l4| , respectively. We can see that the self-energy 



AE from (42) is smoother than in the weak-coupling theory, since the imaginary part of the total self-energy does not 
vanish at the Fermi energy. The mean-field self-energy E — AE is dominant in the vicinity of the Fermi energy, while 
AE, due to the ring diagrams, is responsible for the band tails. 

VIII. CONCLUSIONS 

In this paper we derived a general expansion scheme around mean-field theories for interacting electrons. We used 
a quantum version of the linked cluster expansion modified for the needs of perturbing the mean-field theories under 
consideration. We used the difference of the full and the mean-field Hamiltonian as the perturbation generating the 
LCE. We succeeded in this way to put expansions around weak-coupling (Hartree-Fock) and strong-coupling (Hubbard- 
III) mean-field theories on the same footing. Moreover, the proposed expansion scheme applies to the grand potential 
(free energy) and hence the self-consistent theory has a variational character and is thermodynamically consistent and 
conserving in the Baym sense. The presented application of the LCE around mean fields enables to systematically 
classify, renormalize and sum various classes of diagrams as is usual in the weak-coupling expansion, although the 
Wick theorem and the Gaussian integration are not used. 

We explicitly investigated an approximation using a sum of the so-called ring diagrams leading at the classical 
level to the Bloch-Langer approximation [p!s| . We showed that this approximation at the quantum level, where 
additive internal degrees of freedom represented by Matsubara frequencies are to be involved, is not solvable. Further 
approximations are enforced to obtain the grand potential in closed form. We discussed and compared various 
approximations resulting from simplifications of the complete, renormalized sum of the ring diagrams. Thereby we 
paid special attention to the problem of renormalizations and the self-consistency of the derived approximate solutions. 

At first we studied the influence of dynamical fluctuations due to the U 2 term. When expanded around the 
Hartree-Fock solution, we obtained a modified iterated perturbation theory of ]lO[] . The difference to the IPT was 
manifested by a complex prefactor in the self-energy. This prefactor is a consequence of a partial self-consistency of 
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the approximation and of the fact that the LCE is a diagrammatic expansion for the grand potential, whereas the 
IPT for the self-energy. These prefactors are ubiquitous feature of non-self-consistent or only partially self-consistent 
theories, where the Ward identities and sum rules need not be fulfilled. Although there are no unphysical complex 
prefactors in the IPT for the self-energy, the simple analytic structure of e.g. susceptibilities goes lost. There is no 
set of diagrams generating susceptibilities in the IPT. We found that incompletely self-consistent theories containing 
the U 2 term break down numerically at a critical U c above which no stable numerical solution for complex energies 
was found. This situation did not improve even if we took into account the ring diagrams without any dynamical 
renormalizations. A pole in the two-particle function destroys convergence of the iterations. 

We concluded from the difficulties with nonrenormalized partial summations of the LCE that only dynamically 
renormalized, fully self-consistent theories can provide global, thermodynamically consistent approximations. We 
hence developed a scheme how to introduce full self-consistency into the LCE in a controlled and systematic way. 
We demonstrated the method on an example of the ring diagrams, where we showed that a class of "noncrossing 
diagrams" , when summed, leads to full self-consistency of the approximation. Such a fully self-consistent theory has a 
diagrammatic expansion for any derivative of the free energy. The self-consistent version of the ring diagrams around 
the Hartree-Fock solution was identified with the bubble-chain approximation studied at the level of the self-energy 
in |l9| ]. We studied approximations with the self-consistent ring diagrams for both the weak and strong coupling 
mean-field theories. In the weak-coupling case we found that for intermediate interaction strengths a pole in the 
two-particle function interferes and hinders to reach convergence of iterations. A more advanced approach to handle 
the behavior of the iterations in the vicinity of the pole must be used to decide how the pole is being approached. 
This problem, tightly connected with the metal-insulator transition, will be discussed in a separate publication. The 
ring-diagram approximation around the strong-coupling MFT was found to be stable for all interaction strengths and 
to lead to a metal-insulator transition at about U c ~ 5.5 in the same way as the static MFT of Q. 

The linked cluster expansion with dynamical renormalizations introduced in this paper is not only suitable to derive 
conserving approximations with dynamical fluctuations beyond static mean-field theories, but it may prove efficient 
for expansions around other nontrivial solutions such as the expansion for the Kondo model proposed earlier by 
Anderson and Yuval (^(| . 

We thank D.Vollhardt for useful discussions. The work was supported in part by the grant No. 202/95/0008 of the 
Grant Agency of the Czech Republic and by the Sonderforschungsbereich 341 of the Deutsche Forschungsgemeinschaft. 



APPENDIX A: ELECTRON-HOLE SYMMETRY 



We now prove that the Green function G(z) and consequently the self-energy transform under the electron-hole 
transformation defined as 
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(Al) 



according to formulas ( |3£ 

We use the results of the static mean-field theory |l4|] and the fact that the ring diagrams do not change the 
equations for the mean-field parameters n s and E to obtain the transformation rules 
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(A2) 



Note that here we have to replace the interaction strength U from (T^] with U/2 in the transformation formulas, since 
we use the convention X a = 1/2, X at — instead of X a = 1, X a t = — L 

To prove the electron- hole symmetry of the ring-diagram approximation we assume the validity of (33) and derive 
transformati ons f or the functions G^ (z) and G^ u ' (z) and the integrals (z) and 1^ (z) and then insert them into 
the r.h.s. of (29a) to show that the l eft a s well as the right hand side transform identically. 

Using the definitions in eqs. (|22[), (|26a|), and (J33[) we easily obtain 



G<°>(z)--G< t/ )(-z) 



G^ u \z)^-G^{-z) 



(A3) 



Further we must find transformation rules for the integrals I^(z). If we extend the integration in these integrals 
onto the whole real axis we obtain 



E 



/oo 
dx\C a {z-x) G^{x + tQ + ) -G^(x-i0+) 
■00 L L 

(z + x) [CJ-x + i0+) - C a (-x - iQ+)] } = 0, 



(A4) 
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which follows from the analyticity of the functions G ( -°- ) (z) and C a (z) in the upper and lower complex half-planes. This 
enables to deform conveniently the integration contour into the complex plane and split it into segments containing 
no singularities. 



Using now (A4) and the transformation rules (A3) we obtain 

/(<>)(*) --/OOf-*), iW)( z ) 



When we insert (A3), (A5) into the r.h.s. of (29a) we obtain 

G{z) -» -G(-z) 



-!<«»(-*) 



(A5) 



(A6) 



for both sides of (29a) and hence the equations of motion obey the electron-hole symmetry. 
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Figure Captions 



Fig.l DOS at weak coupling calculated from the LCE with [/ 2 -corrections to the Hartree-Fock MFT at U =1.4 
(broken line) and U = 1.7 (solid line). We used the semi-elliptic DOS with the bandwidth 2D = 4. 

Fig. 2 DOS at U = 1.9 in the same theory as in Fig.l. For larger [/'s the iterations do not converge any longer. 

Fig. 3 DOS calculated with the iterated perturbation theory (IPT) within the domain of convergence of iterations, 
U = 1.7,2.7,4.7 depicted as dotted, broken and solid lines, respectively. 

Fig. 4a Real part of the propagator G (0) within IPT for U = 1.7,3.7,4.7 as in Fig.3. 

Fig. 4b Imaginary part of within IPT for the same values of U as in Fig. 4a. 

Fig. 5 The location oj m i n of the minima in ReG^°^(w) (broken line) and ImG'°'(w) (solid line). 

Fig. 6 DOS for the IPT extended of the ring diagrams at U = 1.87. 

Fig. 7 Simplest diagrams contributing to the expansion of the exponential in (^Oj). The hatched bubbles represent the 
function C, the wavy lines derivatives w.r.t. generating sources. 

Fig. 8 DOS for the fully self-consistent approximation with the ring diagrams around the Hartree-Fock MFT at 
U = 1.5 (broken line) and U = 3.5 (solid line). 

Fig. 9 DOS for the fully self-consistent approximation with the ring diagrams around the strong-coupling MFT at 
U = 1.5 (broken line), U = 3.5 (dotted line), and U = 5 (solid line). 

Fig. 10a Contribution to the imaginary part of the self-energy AE due to the ring diagrams at U — 3.5 for the 
strong-coupling (broken line) and the Hartree-Fock (solid line) MFT. 

Fig. 10b Contribution to the imaginary part of the self-energy E — AE due the Hubbard-III term in the static 
strong-coupling MFT (broken line) and the dynamic MFT with ring diagrams (solid line) at U = 3.5. 
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Fig.7 




Fig. 10a 
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